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SECTION  1 
INTRODUCTION 


Plasmas,  in  undisturbed  or  disturbed  ionospheres  and  magnetospheres, 
can  have  an  enormous  range  of  parametric  values.  Accordingly,  to  understand 
the  development  of  structures  (e.g.,  striations)  under  all  possible 
conditions,  it  is  necessary  to  develop  a  formulation  which  includes  a 
sufficient  amount  of  physics  to  permit  general  applicability.  Laboratory 
simulations  of  HANE  environments  can  contribute  substantially  to  the 
development  of  a  complete  phenomenological  understanding  of  structure 
evolution  in  high-altitude  nuclear  plasmas  as  physical  processes  occur  in  a 
complete  and  natural  way  in  the  laboratory.  Indeed,  laboratory  simulations 
could  provide  crucial  and  objective  guidance  in  determining  those  physical 
effects  which  should  be  carefully  considered  in  analytical  and  numerical 
models.  By  contrast,  subjective  decisions  are  often  made  in  analytical  and 
numerical  studies  in  which  plasma  models  are  often  predicated  on  simplicity  or 
computational  convenience. 

As  an  example  of  how  laboratory  simulations  may  be  used  to  improve 
analytical  and  numerical  models,  consider  the  possibility  that  laboratory 
simulations  might  demonstrate  that  kinetic  effects  (e.g.,  anomolous  diffusion, 
drift  waves,  etc.)  crucially  affect  structure  evolution  in  some  HANE  environ¬ 
ments.  For  such  circumstances,  the  laboratory  measurements  could  suggest 
effective  kinetic  transport  coefficients  for  use  in  fluid  simulations  which 
cannot  consider  kinetic  processes  by  themselves. 

Furthermore,  the  laboratory  simulations  might  demonstrate  that  finite 
temperature  effects  (e.g.,  temperature  gradients  and  thermal  conductivity)  or 
the  detailed  dynamics  of  neutrals  (e.g.,  neutral  viscosity  or  the  interaction 
of  ions  and  electrons  with  neutrals)  are  important  in  some  HANE  environments 
and  should  be  included  in  fluid  simulations  of  later-time,  high-altitude, 
nuclear,  plasma  structures.  By  contrast,  present  fluid  simulations  of  later- 
time  plasma  structures  generally  do  not  include  the  contribution  from  plasma- 
temperature  equations  or  equations  describing  neutral  dynamics.  Hence, 
laboratory  simulations  of  HANE  structure  can  have  a  very  practical  and  crucial 
impact  in  the  objective  determination  of  the  critical  subset,  within  the 
complete  set  of  plasma  processes,  which  should  be  carefully  considered  in 
analytical  and  numerical  models  of  the  HANE  environment. 
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Through  their  impact  on  the  development  of  appropriate  analytical  or 
numerical  models  and  through  the  data  they  provide,  laboratory  simulations  of 
HANE  plasma  structures  could  provide  the  critical,  but  presently  unavailable, 
experimental  information  for  the  development  of  accurate  microstructure 
algorithms  applicable  to  nuclear-effects  simulations  like  SCENARIO. 

The  general  objective  of  the  research  was  to  determine  the  feasibility 
of  using  a  laboratory  device  to  simulate  later-time  high-altitude,  nuclear- 
structured  plasmas.  In  this  work,  we  developed  scaling  laws  between  HANE 
environments  and  the  laboratory,  and  we  determined  optimum  experimental 
configurations  for  the  laboratory  simulations.  We  also  developed  adequate 
diagnostic  techniques  to  evaluate  the  simulations. 

Scaling  laws  are  derived  for  relating  the  evolution  of  plasma  struc¬ 
tures  in  the  later-time  HANE  environment  and  the  large  plasma  chamber  at 
JAYCOR.  Because  the  evolution  of  structure  in  later-time  environments  is 
expected  to  be  well  describable  as  an  electrostatic  and  magnetostatic 
phenomenon,  the  proper  scaling  of  the  l'on-gyroradius  is  argued  to  be  extremely 
significant.*  Indeed,  theoretical  analysis  and  more  recent  nonlinear, 
multidimensional  computations  clearly  demonstrate  the  great  significance  of 
finite  ion-gyroradius  effects  on  striation  "freezing".2’3  To  assure  the 
proper  scaling  of  the  ion-gyroradius  for  reasonable  magnetic-field  strengths 
in  the  laboratory,  we  suggest  that  a  fill  gas  with  very  low  atomic  mass  (e.g., 
hydrogen)  may  be  advantageously  used  in  individual  laboratory  simulations. 

The  initial  plasma  configuration  in  the  laboratory  device  and  possible 
laboratory  studies  are  described.  For  the  cylindrical  plasma  chamber  at 
JAYCOR  a  cylindrical  density  profile,  with  the  electron  density  greater  at 
smaller  radii  than  at  larger  radii,  is  desirable  and  experimentally 
reasonable.  The  externally  applied  magnetic  field  is  directed  axially,  and 
the  externally-applied  neutral  wind  is  directed  across  the  cylindrical  plasma. 
With  this  plasma  configuration,  laboratory  simulations  can  be  undertaken  for 
elucidating  the  role  of  parallel  plasma  currents  and  perpendicular  plasma 
processes  on  the  formation  of  striations.  Experiments  can  also  be  conducted 
for  determining  the  evolutionary  characteri sties  of  single  striations  and 
aggregates  of  striations. 

The  plasma  diagnostics,  for  the  laboratory  simulation  of  structured 
plasmas  in  later-time  HANE  plasmas,  are  described  as  a  variety  of  miniaturized 
particle  collecting  Drobes  to  measure  plasma  density,  electron  and  ion 
temperatures ,  and  electron  density  fluctuations. 
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Experiments  have  been  successfully  conducted  to  demonstrate  the 
capability  of  two-probe  correlation  techniques  to  determine  power-spectral 
densities  in  plasmas  with  density  fluctuations. 

The  well-established  technology  of  Laval  nozzles  is  described  as  the 
means  for  obtaining  the  required  neutral  winds  in  the  laboratory. 

This  report  is  divided  into  seven  sections.  Section  2  derives  scaling 
laws  relating  the  high-altitude  nuclear  environment  and  the  laboratory  tank  at 
JAYCOR.  Section  3  describes  the  initial  plasma  configuration  in  the 
laboratory  device  and  possible  laboratory  studies.  Section  4  describes  the 
diagnostics,  and  Section  5  describes  how  two-probe  correlation  measurements 
can  be  used  to  determine  power-spectral  densities.  Section  6  describes  how  a 
Laval  nozzle  can  be  used  to  obtain  the  required  neutral  winds.  Section  7  is  a 
summary  and  gives  concluding  remarks. 


SECTION  2 
SCALING  LAWS 


Scaling  laws  are  derived  for  relating  the  evolution  of  plasma 
structures  in  the  later-time  HANE  environment  and  the  large  plasma  chamber  at 
JAYCOR.  Consideration  is  given  to  the  possibility  that  the  laboratory 
simulations  may  benefit  from  using  matter  different  from  that  present  in  the 
HANE  environment.  Throughout,  the  symbol  —  designates  the  ratio  of  laboratory 
to  high-altitude  nuclear  quantities. 

In  the  ionosphere  the  collisional  drag  between  a  neutral  wind  and  ions 
is  a  significant  driver  in  striation  formation;  consequently,  it  is  desirable 
to  preserve  the  ratio  between  mean  free  path  and  characteristic  distance: 

r  =  A  .  (1) 

a  L 

In  Eq.  (1)  the  following  symbols  are  used:  p  (mass  density),  A  (atomic  mass), 
a  (binary  collision  cross  section),  and  L  (characteristic  distance).  Ion  and 
neutral  velocities  are  allowed  to  scale  like  the  ion-thermal  speed: 


with  V  (velocity)  and  T  (temperature) ;  therefore,  a  reasonable  representati on 
for  the  scaling  of  time  (t )  is 


From  Eqs.  (1)  and  (3)  it  follows  that  the  binary  collision  frequencies  (v) 
scale  properly  with  time: 


Two  different  ways  are  considered  for  scaling  the  magnetic  field  (B). 
First,  if  the  Alfven  speed  is  preserved,  then 


f  ,  T  s  T  °-5  .  JL 


x  \  0.5 


Second,  if  the  ion-gyroradius  is  preserved,  then 


i.IfiTIs£iC  ,  (6) 

G  QL  Q r 

with  IJ  being  the  scaling  factor  for  charge  state.  The  two  different  scalings 
for  the  magnetic  field  are  completely  consistent  only  if 
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Equation  (7)  is  independent  of  the  scaling  factor  for  atomic  mass,  A.  The 
discrepancy  between  Eqs.  ( 5 ) - ( 6 )  also  occurs  in  scaling  studies  between  the 
early-time  HANE  environment  and  the  laboratory.*’4 

In  the  directions  transverse  to  the  magnetic  field,  striation  evolution 
in  the  later-time  HANE  environment  is  expected  to  be  an  electrostatic  and 
magnetostatic  phenomenon.5,6  Consequently,  plasma  dynamics  associated  with 
the  propagation  of  Alfven  waves  across  the  magnetic  field  should  have  a 
corresponding  negligible  importance  on  the  striation  evolution.  This 
contrasts  with  the  much  more  rapid  early-time  evolution  of  the  HANE 
environment  when  the  geomagnetic  field  is  strongly  modified,  and 
electromagnetic  effects  (e.g.,  Alfven  waves)  are  important.4 

Provided  that  all  electric  fields  have  a  predominantly  electrostatic 
character,  Eq.  (6)  is  easily  shown  to  demonstrate  proper  scaling  between 
Pedersen  and  parallel  conductivities;  therefore,  there  is  proper  scaling 
between  perpendicular  and  parallel  electric  fields.  However,  with  the  use  of 
Eq.  (6),  the  propagation  speed  for  Alfven  waves  along  the  magnetic  field  is 
much  more  rapid  in  the  laboratory  than  if  Eq.  (5)  is  used  instead  to  determine 
the  magnetic-field  strength.  This  means  that  parallel  electric  fields 


correspondingly  tend  to  be  more  electrostatic  in  the  laboratory  than  in  the 
high-altitude  nuclear  environment.6,6  In  this  regard,  theoretical  studies  at 
JAYCOR  suggest  that  the  formation  of  striations  occurs  for  scale  sizes 
transverse  to  the  magnetic  field  which  are  sufficiently  small  so  that  parallel 
currents,  whether  inductive  or  electrostatic  in  character,  are  of  secondary 
importance. ^ »6 

In  contrast  to  the  Alfven  speed,  the  finite  size  of  ion  gyroradii  has 
been  demonstrated  to  have  a  very  strong  effect  on  the  linear  evolution  of 
collisional  Rayleigh-Taylor  instabilities;  ’  moreover,  the  proper  scaling  of 
diffusion  transverse  to  the  magnetic  field  requires  that  the  ion  gyroradii,  in 
addition  to  collision  frequencies,  scale  correctly.  The  conclusion  is  that 
Eq.  (6)  is  much  better  suited  for  scaling  the  magnetic  field  in  the  later-time 
HANE  environment  than  Eq.  (5). 

Figure  1  is  a  summary  of  the  basic  scaling  relations  between  striation 
phenomena  in  the  later-time  HANE  environment  and  the  laboratory. 

Figures  2a  and  2b  are  plots  of  p,  A/T,  and  B„  versus  L  for  T  =  a  =  Q  = 
_  b 

1.  The  scaling  factor  is  included  in  the  figures  for  completeness.  The 
difference  between  Figures  2a  and  2b  is  that  the  gas  used  in  the  laboratory  is 
assumed  to  be  air  (A  =  1)  and  hydrogen  (A  =  0.05),  respecti vely.  For  both 
Figures  2a  and  2b,  it  is  uniformly  true  that  BQ  >  B^;  consequently,  the 
constraint  that  the  ion-gyroradius  scale  properly  results  in  larger  magnetic 
fields  in  the  laboratory  than  would  be  necessary  if  the  scaling  of  Alfven 
waves  happened  to  be  appropriate.  The  use  of  a  gas  with  reduced  atomic  mass 
in  the  laboratory  (i.e.,  A  <  1)  does  have  the  significant  advantage  of 
permitting  reduced  magnetic-field  strengths  in  the  laboratory  consistent  with 
the  scaling  laws. 

As  specific  quantitative  examples,  consider  a  characteristic  scale  size 

of  ~  104  cm  and  ~  1  cm  for  striations  in  the  HANE  environment  and  the 

-  -4 

laboratory,  respectively;  hence,  L  =  10  .  Figure  2a  indicates  that  the  mass 

density,  evolution  rate,  and  magnetic-field  strength  are  104  times  larger  in 
the  laboratory  than  in  the  HANE  environment.  If  the  magnetic-field  strength 
is  0.3  G  in  the  HANE  environment,  then  the  scaling  law  for  ion  gyroradius 
suggests  that  the  magnetic  field  be  3  x  10^  G  in  the  laboratory.  Analogously 
from  Figure  2b,  L  =  104  implies  that  the  mass  density  is  5  x  10^  times  larger 
in  the  laboratory  than  in  the  HANE  environment.  In  the  laboratory,  the 
evolution  rate  and  the  magnetic-field  strength  are  4.5  x  104  and  2.2  x  10^ 
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mass  density: 


A 


velocity: 


time: 


magnetic  field: 


B  ■  <!-*>' - 
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with 


A  *  scaling  factor  for  atomic  mass 
o’  *  scaling  factor  for  collisional  cross  section 
r  *  scaling  factor  for  characteristic  distance 
T  *  scaling  factor  for  temperature 
"Q  *  scaling  factor  for  charge  state 


times  larger  than  in  the  HANE  environment,  respectively.  For  a  magnetic-field 
strength  of  0.3  G  in  the  HANE  environment,  a  magnetic-field  strength  of  6.6  x 
10^  G  is  implied  in  the  laboratory. 

Since  magnetic  field  strengths  <  3  x  10^  G  should  be  attainable  with 
the  JAYCOR  plasma  chamber,  air  can  be  a  realistic  and  reasonable  material  for 
at  least  some  of  the  laboratory  experiments.  In  addition,  the  use  of  a  low 
atomic  mass  material  like  hydrogen  should  permit  parametric  studies  over  an 
especially  large  range  of  L,  the  scaling  factor  for  characteristic  distance. 


SECTION  3 

INITIAL  PLASMA  CONFIGURATIONS  AND  LABORATORY  STUDIES 

The  initial  plasma  configuration  shown  in  Figure  3  appears  reasonable 
for  the  laboratory  simulation  of  structured  plasmas  in  the  later-time  HANE 
environment.  A  cylindrical  electron-density  profile,  n(r),  is  created  in  the 
JAYCOR  plasma  chamber  tank  with  the  electron  density  being  greater  at  smaller 
radii  than  at  larger  radii.  This  density  profile  can  be  created  by  injection 
and  subsequent  rapid  ionization  of  neutrals.  The  density  profile  does  not 
need  to  be  centered  on  the  axis  of  the  plasma  chamber.  The  axial,  z-directed, 
magnetic  field  with  strength  B  confines  electrons  and  ions.  With  a  neutral 
wind,  vp,  directed  across  the  plasma,  the  density  gradient  and  the  neutral 
wind  are  antiparallel  on  the  backside  of  the  cylindrical  plasma;  consequently, 
the  plasma  becomes  unstable  to  the  Rayleigh-Taylor  modes  which  are  thought  to 
be  primarily  responsible  for  the  formation  of  striated  plasmas  in  the  HANE 
envi ronment. 

Even  with  the  constraints  imposed  by  the  scaling  laws  described  in 
Section  2,  the  basic  plasma  configuration  has  substantial  flexibility  to 
permit  various  types  of  relevant  experiments.  In  the  following  paragraphs, 
the  different  types  of  experiments  are  synopsized  through  their  division  into 
three  categories:  (1)  parallel  plasma  currents  and  the  formation  of 

striations,  (2)  perpendicular  plasma  processes  and  the  formation  of 
striations,  and  (3)  evolution  of  single  and  aggregates  of  striations. 

3.1  PARALLEL  PLASMA  CURRENTS  AND  THE  FORMATION  OF  STRIATIONS. 

References  6-8  provide  clear  quantitative  demonstrations  of  how 
structures  with  large  scale  sizes  across  the  magnetic  field  map  greater 
distances  along  the  magnetic  field  than  structures  with  small  transverse  scale 
sizes.  The  key  physical  process  is  the  short-ci rcuiting  of  differences  in  the 
electrostatic  potential  between  adjacent  magnetic-field  lines  by  the  combina¬ 
tion  of  finite  parallel  and  finite  perpendicular  currents.  The  observation  of 
sharp  edges  on  ionospheric  clouds  transverse  to  the  geomagnetic  field9"^  and 
additional  theoretical  evidence^"®  suggests  that  whereas  the  evolution  of 
large  plasma  clouds  may  strongly  depend  on  parallel  currents,  the  steep 
density  gradients  and  the  initial  formation  of  Rayleigh-Taylor  instabilities 
upon  the  steep  density  gradients  are  only  weakly  dependent  upon  parallel 
currents.  Indeed,  the  suggestion  has  been  made  that  an  upper  bound  for  the 


Figure  3.  Initial  plasma  configuration  for  the  laboratory  simulation  of 
structured  plasmas  In  the  later-time,  high-altitude,  nuclear 
environment. 
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scale  size  of  striations  transverse  to  the  magnetic  field  is  determined  by  the 

constraint  that  the  contributions  from  parallel  electron  currents  be 

o 

negligible  during  the  initial  Rayleigh-Taylor  evolution. 

The  use  of  end  plates  with  varying  degrees  of  electrical  insulation 
could  be  used  to  experimentally  determine  how  parallel  plasma  currents  impact 
the  formation  of  sharp  transverse  density  gradients  and  the  subsequent 
formation  of  striations.  With  insulating  end  plates,  parallel  currents  should 
be  minimized;  therefore,  the  transverse  density  gradients  are  predicted  to 
have  their  largest  possible  value,  and  the  maximum  permissible  transverse 
scale  sizes  for  striations  should  increase. 

At  high-geomagnetic  latitudes,  particle  precipitation  and  coupling 
between  the  ionosphere  and  the  magnetosphere  can  substantially  impact  the 
development  of  plasma  structures.  With  the  capability  for  biasing  the  end 
plates  in  the  laboratory  tank  at  different  voltages  and  the  capability  for 
injecting  electron  beams  along  the  externally  applied  magnetic  field, 
laboratory  simulations  of  high-latitude  plasma  structures  are  possible.  By 
biasing  the  end  plates,  experimental  studies  of  how  applied  plasma  currents 
affect  the  development  of  plasma  structures  can  be  conducted.  Moreover,  by 
injecting  electron  beams  along  the  applied  magnetic  field,  the  effects  of 
localized  particle  precipitation  could  be  studied  as  a  seed  or  as  a  direct 
source  for  plasma  structures. 

If  the  driving  force  transverse  to  the  magnetic  field  (e.g.,  ion- 
neutral  drag  with  a  neutral  wind)  has  finite  extent  along  the  magnetic-field 
direction  as  expected  in  the  later-time  HANE  environment,  the  self-consistent 
current  patterns  are,  in  general,  expected  to  have  some  finite  components 
along  the  direction  of  the  magnetic  field.  In  the  laboratory,  the  effects  of 
finite  parallel  extent  for  driving  forces  can  be  simulated  very  simply  by 
injecting  the  neutral  wind  across  the  magnetic  field  over  a  specified  parallel 
distance.  This  would  permit  study  of  how  the  finite  parallel  extent  of 
driving  forces  impacts  the  development  of  sharp  transverse  density  gradients 
and  the  subsequent  initial  formation  of  striations. 

3.2  PERPENDICULAR  PLASMA  PROCESSES  AND  THE  FORMATION  OF  STRIATIONS. 


Various  theoretical  and  numerical  studies  have  examined  how  plasma 
processes  acting  in  the  directions  across  the  magnetic  field  affect  the 
characteristics  of  striations,  especially  their  transverse  scale  sizes. 


Diffusion5-11  and  ion  viscosity^  are  some  of  the  more  important  plasma 
processes  which  have  been  considered. 

As  an  example  of  how  the  laboratory  simulations  can  be  used  to 
determine  the  dominant  processes  affecting  the  evolution  of  striations  in  the 
directions  transverse  to  the  magnetic  field,  consider  the  suggestion  of 
Reference  2  that  the  lowest-order  ion-gyroradius  contributions  to  the  ion- 
viscosity  tensor  are  the  major  contributors  to  the  observed  scale  sizes  of 
striations  transverse  to  the  magnetic  field.  In  particular,  the  theory  of 
Reference  2  proposes  the  following  expression  for  the  tranverse  wave  number, 
kt,  characteristic  of  "frozen"  striations: 


1  /V-in\°*i 

k  — 

1  pi  \fli  / 


In  Eq.  (8),  the  following  symbols  are  used:  (ion-gyroradius),  v-j n  (ion- 

neutral  collision  frequency),  and  (ion  gyrofrequency).  Figure  1  in  Section 
2  provides  the  relationship  between  parameters  in  the  later-time  HANE 
environment  and  the  laboratory.  The  application  of  the  scaling  laws  to  Eq. 
(8)  indicates  that  kt  scales  properly  as  inverse  distance;  therefore,  Eq.  (8) 
not  only  gives  a  prediction  for  transverse  striation  scale  sizes  in  the 
laboratory  but  it  also  predicts,  in  conjunction  with  the  scaling  laws, 
corresponding  transverse  striation  scale  sizes  appropriate  to  the 
ionosphere.  The  observation  that  Eq  (8)  is  satisfied  in  the  laboratory  would 
provide  strong  experimental  evidence  .or  its  applicability  in  the  later-time 
HANE  environment  within  the  ionosphere. 


3.3  EVOLUTION  OF  SINGLE  AND  MULTIPLE  STRIATIONS. 

Figure  4  illustrates  the  qualitative  evolution  of  striations  as  deduced 
from  the  observation  of  ionospheric  barium  clouds.^  Similar  evolution  may  be 
expected  for  HANE  plasmas  within  the  ionosphere.  In  particular,  the  initial 
evolution  of  linear  Raylei gh-Taylor  instabilities  is  soon  followed  by 
nonlinear  development  into  sheets  of  ionization  aligned  along  the  direction  of 
the  magnetic  field.  The  sheets  themselves  can  evolve  in  a  complex  way  as  they 
may  bifurcate,  or  they  may  form  rods  as  they  pinch  away  from  the  bulk  of  the 
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Figure  4.  Qualitative  evolution  of  striatlons  as 
observation  of  Ionospheric  barium  clouds.9 


deduced 


plasma  cloud.  As  time  progresses,  the  number  of  rods  increases.  The  rods 
eventually  dissipate  as  they  ablate  into  very  thin  sheets. 

A  successful  laboratory  simulation  of  striation  evolution  should  also 
be  expected  to  have  the  qualitative  features  illustrated  in  Figure  4.  A  major 
advantage  of  the  laboratory  simulation  would  be  that  appropriately  placed 
wires  could  be  used  to  affect  the  development  of  individual  striations.  For 
example,  with  the  capability  of  being  able  to  eliminate  all  striations  but  one 
with  wire  barriers,  experiments  could  be  conducted  to  study  how  the  presence 
of  neighboring  striations  affects  the  development  of  individual  striations. 
Also,  if  wires  could  be  used  to  destroy  the  plasma  bridge  between  the  tips  of 
striations  and  the  bulk  plasma,  then  experiments  should  be  possible  clarifying 
the  significance  upon  striated  plasmas  of  attachment  or  disattachment  to  bulk 
plasma. 

For  nuclear-effects  simulations  like  SCENARIO, 12-13  characteristic 
scale  sizes  for  striations  transverse  to  the  magnetic  field  are  much  smaller 
than  the  distance  corresponding  to  the  spacing  between  grid  points.  The 
striations  are  microscopic  features  for  the  simulations,  and  effects 
attributable  to  the  presence  of  striations  must  be  treated  in  a  statistical 
manner.  The  striation  convection  algorithm  in  SCENARIO  calculates  the 
convection  of  striations  in  a  statistical  way  throughout  the  later-time,  high- 
altitude,  nuclear  environment.12  An  implicit  but  unproven  assumption  in  the 
derivation  of  the  algorithm  is  that  the  motion  of  individual  striations  is 
determined  by  field-line-integrated,  transverse,  current-densi ty  fields  over 
scale  sizes  much  larger  than  the  transverse  scale  size  of  individual 
striations;  therefore,  the  proximity  of  a  single  striation  to  its  nearest 
neighbors  is  assumed  to  have  no  special  significance. 

In  the  laboratory  simulation,  the  neutral  wind  can  be  permitted  to  have 
shear  in  the  directions  transverse  to  the  magnetic  field;  therefore,  the 
motion  of  fully-formed  striations  can  be  studied  to  determine  whether  they 
move  in  a  manner  consistent  ith  the  direction  of  the  neutral  wind  and  large- 
scale,  self-consistent,  electric-field  forces.  Furthermore,  if  wire  barriers 
can  be  used  to  eliminate  some  striations,  it  should  be  possible  to  determine 
the  degree  that  proximity  to  neighboring  striations  affects  striation 
convection. 


SECTION  4 

DIAGNOSTICS  FOR  THE  LABORATORY  SIMULATIONS 


Plasma  diagnostics,  for  the  laboratory  simulation  of  structured  plasmas 

in  later-time  HANE  plasmas,  consist  of  a  variety  of  miniaturized  particle 

collecting  probes  to  determine  plasma  density,  electron  and  ion  temperature, 

and  density  fluctuations.  Cylindrical  Langmuir  probes  are  used  to  monitor 

plasma  density  and  T  ,  and  they  can  also  be  the  primary  method  for  measuring 

i  ,e 

the  temporal  and  spatial  development  of  plasma  structures,  e.g.,  striations 

during  the  neutral  wind  injection.  Single-sided  planar  disc  probes  will  be 

used  to  monitor  T  .  Two  of  the  disc  probes  can  be  mounted  back-to-back  and 
i  ,e 

monitored  differentially  to  determine  if  any  net  current  flow  exists  along 
Bq.  To  measure  the  ion  temperature,  T^,  a  gridded,  biased  Faraday  cup  (energy 
analyzer)  can  be  constructed.  A  schematic  drawing  of  such  an  ion  analyzer  is 
given  in  Figure  5.  Principles  of  operation  and  specific  design  details  may, 
for  all  three  of  the  above  mentioned  probes,  be  found  throughout  the  plasma 
physics  literature.14"19 

One  difficulty  involved  in  measuring  the  temperature,  T,  for  both  ions 
and  electrons  is  the  interpretation  of  probe  currents  in  a  high  magnetic 
field.  For  BQ  s  3  kG,  the  Larmor  radius  for  the  electrons  will  be  rg  <  10"3 
cm,  which  is  much  smaller  than  the  anticipated  probe  dimensions  of  0.1  cm  x 
10"2  cm  (radius).  The  electron  gyrof requency ,  uce/2n  =  10  GHz,  may  also  be 
greater  than  the  electron  plasma  frequency  u>pe/2n  =  1  GHz  at  ne  =  1010  cm'3. 
Thus,  the  electrons  can  be  strongly  magnetized;  therefore,  to  inter¬ 
pret  T  from  the  cylindrical  probe  I-V  characteristic  must  involve  scaling 
the  temperature  obtained  from  a  more  exact,  albeit  more  complicated  metnod, 

such  as  using  test  ion  acoustic  waves.  By  monitoring  the  wave  length  at  a 
given  driving  frequency  greater  than  the  ion-cyclotron  frequency,  the  ion 

sound  speed  can  be  determined  as  Cs  =  fx .  Then,  using  the  relation  for  these 

O 

waves,  =  kTe/M^ ,  and  the  known  ion  mass,  a  temperature  may  be  inferred. 
Temperatures  measured  and  compared  using  these  two  techniaues  have  bee-'  show'' 
to  differ  by  not  more  than  10-15%  in  the  present  configuration  wt*  i*e 

magnetic-field  strength  of  50  gauss. 

Absolute  density  measurements  using  Langmuir  probes  are  subject  to 

large  errors  in  a  magnetic  field.  Thus,  only  relative  density  and  fluctuation 
measurements  can  be  made  using  these  probes.  "he  absolute  sens' ty  can  ne 
measured  by  observing  the  cone-angle  of  rf  nesjnan^  surfaces  excited  witn  3 


small  rf  probe.  From  the  cone  angle  Oc,  the  electron  plasma  frequency  upe 
is  computed  according  to 


iu2  =  (u>2 

pe  '  ce 


2  .  2 
“ces,n  ec 


knowing  wCP  and  the  rf  driving  frequency  u>. 

To  determine  u>ce,  or  essentially  B,  diamagnetic  loop  probes  can  be  used 
to  monitor  the  time  varying  magnetic  field  3B/3t.  It  may  also  be  possible  to 
use  a  heat  resistant  Hall  probe  with  moderate  time  response  (10  kHz)  to 
measure  the  field  directly. 


SECTION  5 

TWO-PROBE  CORRELATION  STUDIES 


The  measurement  of  power-spectral  densities  for  the  scaled  laboratory 
simulations  should  provide  a  significant  basis  for  determining  the  effects  of 
actual  HANE  structures  on  communications  and  tracking;  therefore,  we  have 
undertaken  actual  experiments  to  demonstrate  our  capability  to  measure 
frequency  and  wave  number  spectrums. 

For  these  experiments,  plasma  turbulence  was  generated  by  an  axial 

electron  current  (parallel  to  B0)  which  thereby  excited  a  broad  spectrum  of 
ion  acoustic  and  electron  cyclotron  waves  in  the  background  steady-state 
plasma.  The  electron  drift  was  driven  by  a  positive  voltage  pulse  applied  to 
a  15  cm  diameter  anode  plate  at  one  end  of  the  plasma  chamber.  The  pulse 
width  was  varied  from  1  to  10  ms  but  was  usually  set  at  4  ms.  This  method  of 
generating  the  turbulent  wave  spectrum  was  used  since  it  was  found  that  the 
turbulence  inherent  in  the  plasma-density  gradients  at  low  magnetic  field 
values  (Bq  >  10  gauss)  was  such  that  the  correlation  length  rc  was  much 

smaller  than  a  typical  wavelength.  Thus,  a  more  controlled  means  of  producing 

the  waves  was  used,  yielding  turbulent  waves  that  were  correlated  to  distances 
comparable  to  a  wavelength. 

The  two  probes  consisted  of  2  mm  diameter  tantalum  discs  (A  =  2ur2  = 
6.3  x  10"2  cm^)  mounted  axially  and  radially  with  respect  to  the  main 

turbulent  plasma  column  which  had  a  diameter  roughly  defined  by  the  anode 
plate  (15  cm)  and  an  axial  extent  of  40-50  cm  in  the  region  of  uniform 
magnetic  field.  Both  probes  were  mounted  on  precision  drive  mechanisms  which 
could  move  in  0.5  mm  ±  0.05  mm  steps.  The  measured  turbulent  spectrum  is 
shown  in  Figure  6,  as  6n/n  versus  frequency  f.  This  was  swept  over  the 
frequency  range  of  interest  (10  kHz-2  MHz). 

A  schematic  of  the  circuit  used  to  make  the  wavelength  measurements  of 
the  turbulent  waves  is  shown  in  Figure  7.  The  probe  was  essentially  monitor¬ 
ing  fluctuations  in  the  electron  density  associated  with  the  ion  acoustic 
turbulence.  This  was  done  by  biasing  the  probe  just  above  the  plasma 
potential  and  measuring  the  fluctuating  electron  current  drawn  to  the  probe  as 
a  voltage  drop  across  the  1  kn  load  resistor.  This  voltage  signal  was  then  ac 
coupled  into  a  preamplifier  (gain  of  100)  and  then  fed  to  narrow  (Wu>  -  1/35) 
bandpass  filters.  A  schematic  for  these  "active"  filters  which  were  con¬ 
structed  and  calibrated  at  JAYCQR  is  given  in  Figure  8.  A  plot  of  the 
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Figure  6 


Measured  frequency  spectrum  of  ion  acoustic  turbulence 
Hg  plasma,  n  =  2  x  109  cm-3,  1^  =  0.20  A. 


Schematic  of  the  two-probe  autocorrelation  configuration. 


frequency  response  of  a  filter  tuned  to  93  kHz  is  shown  in  Figure  9,  measured 
using  a  swept  frequency  network  analyzer.  These  filters  boosted  the  signal 
amplitude  at  the  pass-band  frequency  an  additional  factor  of  10  and  then  a 
unity  gain  buffer  amplifier  was  used  to  drive  the  50  a  impedance  of  the  mixer 
(signal  multiplier). 

The  mixer  output  contained  the  sum  and  difference  frequencies  of  the 
two  probe  signals,  so  a  low-pass  filter  was  used  to  suppress  the  2 u>  portion  of 
the  signal.  This  low  frequency  signal  (~dc  for  test  waves  and  ~  1  kHz  for 
turbulent  waves)  was  then  averaged  over  many  anode  plate  current  pulses  using 
a  digital  storage  oscilloscope  (Tektronix  7854).  With  one  probe  fixed,  the 
other  probe  was  moved  and  the  average  mixed  signal  value  was  recorded  and 
plotted  versus  probe  position  as  an  interferogram  to  determine  the  wavelength 
of  the  particular  frequency  mode. 

Small  amplitude  test  waves  were  launched  using  the  reference  (fixed) 
probe  as  a  check  on  the  electronics  to  ascertain  what  wavelengths  we  could 
expect  to  see  in  the  turbulent  wave  case.  Figure  10  shows  a  test  wave  inter¬ 
ferogram  taken  at  a  driving  frequency  of  f  =  250  kHz.  There  was  a  fairly  good 
correlation  out  to  az  =  20  mm,  and  from  the  measured  wavelength  x  =  8.5  mm, 
the  sound  speed  was  calculated  as  C$  =  fX  =  2.1  x  105  cm/s.  For  the  helium 

plasma  which  was  used  here,  this  corresponds  to  an  electron  temperature  T  = 

2  e 
0.2  eV  since  Cs  =  kTe/M.j.  This  was  comparable  to  the  temperature  measured 

using  a  swept  Langmuir  probe  where  Tg  s  0. 1-0.3  eV.  Figure  11  shows  a  similar 

curve  taken  at  f  =  100  kHz  yielding  Cs  2  1.9  x  105  cm/s  and  Tg  =  0.14  eV.  All 

test  waves  and  turbulent  wave  measurements  were  done  axially,  along  the 

magnetic  field,  since  the  waves  were  quickly  decorrelated  as  one  moved 

radially  across  the  field.  The  anode  plate  was  located  at  z  =  +  200  mm. 

The  turbulent  wave  correlation  data  showed  much  stronger  damping  than 
the  test  wave  case,  probably  due  to  the  much  smaller  amplitude  of  the  wave  (mV 
versus  V).  An  interferogram  taken  at  100  kHz  is  shown  in  Figure  12.  The 
slight  minimum  at  Az  =  0  was  due  to  a  shadow  effect  of  the  probes  being  so 
close  to  each  other.  The  curve  went  through  zero  between  4  and  5  mm  which  was 
interpreted  as  the  X/4  distance,  and  then  rapidly  decayed  to  zero  (no 
correlation).  For  X  2  18  mm,  this  gave  Cs  =  1.8  x  105  cm/s  and  Tg  =  0.13  eV. 
Additional  scans  were  taken  at  250,  500,  750,  and  1,000  kHz  with  similar  wave¬ 
forms.  The  results  at  1  MHz,  however,  were  inconclusive  since  the 
i nterferogram  decayed  monotoni cal ly  to  zero,  with  no  zero  crossings.  (f  ^  ; 


Test  wave  interferogram  for  helium  plasma 
recorded  at  a  driving  frequency  of  250  kHz 


correlation  data  at  100  kHz 


=500  kHz  H  100  SHOT  AVERAGE 


Interferogram  wave  data  recorded  at  750  kHz. 
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SECTION  6 

NEUTRAL-WIND  SOURCE 


To  undertake  laboratory  simulations  of  plasma-structure  formation  in 
the  ionosphere,  neutral  wind  velocities  ~  1  x  103  cm/s  and  neutral  densities 
~  (10^-10*4)  cm"3  are  suggested  by  the  scaling  laws.  The  well-established 
technology  of  Laval  nozzles  can  provide  the  requisite  neutral  winds  with 
velocities  corresponding  to  ~  0.3-1  eV  particle  energies  and  neutral  densities 
~  10n-1013  cm-3;^  therefore,  the  Laval  nozzles  are  clearly  logical  and 
relevant  for  the  laboratory  simulations. 

A  Laval  nozzle  is  an  air-flow  tube  that  is  used  to  produce  supersonic 
gas  flows.  They  are  used  in  wind  tunnel  experiments  and  they  are  used  to 
supply  the  puffed  gas  on  several  of  the  gas  puffed  z  pinch  machines.  Andersen 
et  al.  used  a  mock  Laval  nozzle  to  simulate  in  the  laboratory  supersonic 
expansion  of  a  solar  wind, 20  and  they  formed  virtual  boundaries  on  their 
nozzle  with  a  magnetic  field.  JAYCOR  currently  retains  the  services  of  Will 
Peschel,  an  engineering  specialist  who  has  designed  and  built  hundreds  of 
Laval  nozzles. 

Laval  nozzles  can  readily  supply  high  particle  densities  but  they  are 
limited  to  high  velocities  (greater  than  1  x  105  cm/s)  at  intermediate 
densities  (1  x  1011  to  1  x  1013)  cm"3.  The  nozzles  cannot  function  below  1  x 
1010  cm-3  as  they  must  operate  in  a  pressure  regime  where  the  mean  free  path 
of  the  gas  molecules  is  much  smaller  than  the  smallest  dimension  of  the 
nozzle.  Gas  is  compressed  into  a  Laval  nozzle  and  reaches  a  Mach  1  velocity 
at  the  minimum  nozzle  cross-section.  The  emerging  gas  goes  through  a  pressure 
discontinuity  at  the  Mach  1  interface  and  expands  with  increasing  velocities 
out  through  the  exit.  Mach  numbers  up  to  10  are  easily  achieved.  However, 
the  gas  must  be  in  a  continuum  state  (col  1  isional )  at  the  nozzle  entrance. 
Large  mean  free  paths  inhibit  uniform  gas  flow  and  acceleration.  To  achieve  1 
x  104  cm/s  at  1  x  1011  cm"3  on  the  downstream  side  means  that  we  would  have  to 
use  supercooled  gases  (lower  speed  of  sound)  and  high  Mach  numbers  in  excess 
of  10  (large  area  ratios  for  continuum  gas  flow).  It  is  possible  to  achieve  1 
x  104  to  2  x  10^  cm/s  at  1  x  lO*4  cm-3. 

Figure  16  shows  a  schematic  of  a  possible  laboratory  neutral  wind 
source.  Moving  from  left  to  right  across  the  figure  we  have  a  cooled  neutral 
gas  reservoir,  a  Varian  piezoelectric  puff  gas  valve,  a  flow  quieting  section, 


Figure  16.  Schematic  of  neutral  wind  assembly. 


a  two-dimensional  Laval  nozzle,  our  plasma  vacuum  chamber,  and  an  extension 
drift  tube. 

The  following  specifications  would  be  used  to  build  a  nozzle  capable  of 
delivering  -lx  10^  to  -  2  x  105  cm/s  neutral  wind  at  10"3  torr,  with  a 
uniform  two-dimensional  cross-section  2  x  10  cm  without  flooding  the  vacuum 
chamber.  The  critical  parts  of  the  nozzle  configuration  consist  of  throat 
size  to  molecular  mean  free  path  and  total  gas  injected  compared  with  chamber 

C  *3 

pressure.  A  gas  moving  at  1  x  103  cm/s  will  take  10"°  s  to  “bounce"  off  the 
walls  of  our  chamber.  The  gas  fills  a  2  cm  x  10  cm  x  100  cm  volume  at  10 
torr.  The  flow  rate  equals  the  pressure  divided  by  time  and  multiplied  by  the 
volume.  The  flow  rate  can  also  be  defined  in  terms  of  the  nozzle. 

_  „  K  x  A*  x  P. 


where  w  is  the  flow  rate,  p  is  the  neutral  gas  pressure,  V  is  the  gas  volume 
in  cm,  t  is  its  on  time,  K  is  a  conversion  constant  =  3.4,  A*  is  the  nozzle 
throat  area  in  cm2,  PQ  is  the_  stagnation  pressure  in  torr,  and  T0  is  the 
stagnation  temperature  (determined  from  nozzle  specification  in  Rankines). 
There  are  no  unknowns.  Nozzle  design  tables  dictate  the  throat  area  and 

pi 

stagnation  values. 

From  nozzle  design  tables,  A/A*  (exit  area/throat  area)  equals  4.35  for 
a  Mach  3  nozzle.  P/PQ  (delivered  pressure/stagnation  pressure)  =  2.7  x 
10“2.  Similarly,  T/TQ  =  0.36.  The  speed  of  sound  is  a  =  k'  /T,  a  =  3  x  10^ 
cm/s  *  (Mach  3).  k'  «  50  for  T  in  Rankines. 

P0  =  10_3/2.7  x  10"2  =  0.04  torr 
TQ  =  a2/49  =  500°R  (40°F) 

A*  =  20/4.35  =  4.6  cm2  .  (11) 


Thus,  we  have  the  reservoir  pressure  and  can  infer  that  the  molecule  mean  free 
path  is  small  compared  to  the  throat  diameter.  The  gas  is  maintained  close  tc 
room  temperature  (desired),  and  the  reservoir  pressure  and  nozzle  dimensions 
are  readily  achieveable. 


How  can  we  verify  the  neutral  wind  output?  That  part  is  easier  to 
solve.  We  can  use  ionization  gauge  elements  to  detect  the  fast  gas  flow  and 
quantify  the  neutral  density.  An  ionization  gauge  works  by  ionizing  a  known 
fraction  of  the  surrounding  neutral  gas  population  and  accelerates  the  just 
formed  ions  into  a  collector.  The  gauges  contain  electron  emitting  filaments 
that  ionize  the  neutrals  with  60  eV  electron  energies.  The  gauges  can  respond 
to  microsecond  changes  in  pressure.  Their  response  is  linear  with  pressure  as 
long  as  the  ion  sheath  that  forms  around  the  sensing  element  is  not  space- 
charge  limited.  They  are  sometimes  called  Residual  Pressure  Analyzers  (RPA). 

A  single  gauge  element  cannot  differentiate  between  directed  gas  flow 
and  net  pressure  increase  because  it  measures  the  product  of  density  times 
velocity.  We  would  use  two  gauges  to  monitor  gas  velocity  via  a  time-of- 
flight  measurement.  We  would  separate  the  two  detectors  by  at  least  2  cm  (at 
1  x  105  cm/s  that  means  a  20  microsecond  transit)  to  insure  nonperturbing 
sheath  current  measurements. 
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SECTION  7 
CONCLUSIONS 


Scaling  laws  have  been  derived  for  relating  the  evolution  of  plasma 
structures  in  the  later-time  HANE  environment  with  the  large  plasma  chamber  at 
JAYCOR.  In  the  derivations  it  is  argued  that  the  proper  scaling  of  the  ion- 
gyroradius  is  extremely  important;  therefore,  to  assure  reasonable  magnetic- 
field  strengths  in  the  laboratory,  we  suggest  that  a  fill  gas  with  very  low 
atomic  mass  (i.e. ,  hydrogen)  may  be  advantageously  used  in  individual 
laboratory  simulations. 

The  initial  plasma  configuration  and  possible  laboratory  studies  are 
described.  For  the  cylindrical  tank  at  JAYCOR  a  cylindrical  density  profile, 
with  the  electron  density  greater  at  smaller  radii  than  at  larger  radii,  is 
desirable  and  experimental ly  reasonable.  The  externally  applied  magnetic 
field  is  directed  axially,  and  the  externally-applied  neutral  wind  is  directed 
across  the  cylindrical  plasma.  This  laboratory  configuration  permits  simula¬ 
tions  for  elucidating  the  role  of  parallel  plasma  currents  and  perpendicular 
plasma  processes  on  the  formation  of  striations.  Experiments  can  also  be 
conducted  for  determining  the  evolutionary  characteristics  of  single 
striations  and  aggregates  of  striations.  The  extensive  use  of  probes  to 
measure  plasma  parameters  and  power-spectral  densities  can  assure  the 
relevance  of  the  laboratory  simulations  to  the  understanding  of  plasma 
structures  in  later-time  HANE  plasmas.  The  well-established  technology  of 
Laval  nozzles  can  be  used  to  drive  a  neutral  wind  across  the  cylindrical 
laboratory  plasma. 
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